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ABSTRACT 

In this paper we model the evolution of PrePlanetary Nebula (PPN) and Planetary 
Nebula (PN) morphologies as a function of nebular age. The aim of the work is to 
understand if shape transitions from one evolutionary phase to the other can be driven 
by changes in the parameters of the mass-loss from the central star. We carry out 
2.5D hydrodynamical simulations of mass-loss at the end stages of stellar evolution 
for intermediate mass stars. Changes in wind velocity, mass-loss rate and mass-loss 
geometry are tracked. We focus on the transition from mass-loss dominated by a short 
duration jet flow (driven during the PPN phase) to mass- loss driven by a spherical fast 
wind (produced by the central star of the PN). Our results show that while jet driven 
nebulae can be expected to be dominated by bipolar morphologies, systems that begin 
with a jet but are followed by a spherical fast wind will evolve into elliptical objects. 
Systems that begin with an aspherical AGB wind evolve into butterfly shaped 
nebula with, or without, a jet phase. In addition, our models show that spherical 
nebulae are highly unlikely to derive from either bipolar PPN or elliptical PN over 
relevant time scales. The morphological transitions seen in our simulations may how- 
ever provide insight into the driving mechanisms of both PPN and PN as point to 
evolutionary changes in the central engine. 

Key v^rords: stars: AGB and post- AGB - (ISM:) planetary nebulae: general - stars: 
winds, outflows - hydrodynamics - radiative transfer - methods: numerical 
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1 INTRODUCTION 

Planetary Nebulae (PN) expand at a few lOkms"^ away 
from evolved intermediate-mass stars (with main sequence 
progenitor masses ~1— 8Mq). These nebulae exhibit a 
rich variety of morphologies with the main categories 
being spherical , ellip tical and bipolar (for a review see 
iBalick fc FrankI [2OO2I ). Spherical PN have been explained 
by time-dependent interacting stellar wind models (ISW; 
Kwok, Purton, & Fitzgerald 1978) that describe the neb- 
ula as the collision of a fast (~1000kms~^) tenuous (be- 
tween about 10~* and 10~^ M0yr~^) stellar wind driven 
by the hot central star of PN (CSPN) with a dusty slow 
(~ lOkms"^) dense (~10~^ M© yr~^) wind, or circumstellar 
envelope, ejected during the earlier AGB phase. The nebu- 
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lae are thus formed when the fast wind catches up with the 
AGB wind, driving a shocked shell through it. 

The formation of bipolar PN and jets has been ad- 
dressed by generalized ISW models (GISM) in which the 
ejected AGB envelope is characterized by either a toroidal 
density distribu tion or an as p heric a l velo c ity field (Kahn 
& West 1985; iBalickl 1 19871: Ilck3 1 19881 : iHuggins I l2007l : 
Frank & Mellema 1994: llckd [l988l : Icke, Preston, & Bab 
ick, 1989; Melle ma, Eu lderink, & Icke, 1991; Icke, Balick, 
& Frank, 1992; . Frank fc Mellema 1994,; Dwarkadas, Gheva- 
lier, & Blondin, 1996). In some of these models, asym- 
metries in AGB envelopes develop as a result of a bi- 
nary stellar s ystem, in which an AGB star interacts with 
a companion (Garcia- Ar redondo fc Frank 2 0041 : lEdgar et al.l 
(2008), or an AGB star and a companion share a com- 
mon envelope evolution (Ilbenlll99ll : lLividll993l : ISokeiJll997l : 
ISandguist et al.l 119981 : ISokeiJ Il998l . and references therein; 
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De Marco et all l2003l : iNordhaus et alj l2007l : iPassv et al] 
2012r ). Note, however, that AGB envelopes generally show 



spheri cal shapes (Bujarrabal fc Alcolea 1991; K ahane & Jura 
1994; IStanek et all Il995l : iGroenewegenI Il996l ). Some AGB 
observations show seeds of asymmetry. 

The original GISW model (see e.g. Balick & Frank 2002) 
assumed a spherical fast wind driven from the CSPN ex- 
panding into an aspherical AGB wind. While these models 
were able to generate hot, low density jet flows via shock fo- 
cusing at the inner shock, they were less successful at recov- 
ering dense colder jets like those seen in YSOs. Since many 
PPN and young PN show evidence for such narrow jets, 
Sahai & Trauger (1998) suggested that coUimated PN flows 
created at, or near, the central engine were the real drivers of 
early PN morphology. Such jet formation has been proposed 
as a natural conse quence of interacting binary evolution 
dLee fc Sahail l200l 2004 ; lAkashi fc So"ke3 boosl : iLee et all 
120091 ). iLee fc Sahail l|2003l ) carried out numerical simulations 



and synthetic observations of a jet (or coUimated fast wind) 
interacting with a spherical AGB wind. They found that 
both the dynamics and emission of PPN, and young PN, 
shells depend on the v elocity and geometry o f the jet. Us- 
ing a similar approach lAkashi fc Sokerl l|2008r ) modeled the 
effect of short-duration jets expanding into a spherical, slow 
AGB wind on the observed shape of evolved PN. The sim- 
ulations of lAkashi fc Sokerl (|2008r ) covered tim escales about 
an ord e r of magnitude l o nger than the ones of iLee fc Sahail 
l|2003l ). lAkashi fc Sokeil (|2008l ) found that the AGB mass- 
loss history also has a profound influence on the nebular 
structure, and that the head of the lobes mov e sufficiently 
fast t o excite some visible emission lines (see I Lee fc Sahail 
I2OO3I for detailed results of optical forbidden line emission 
which is excited at the heads of the lobes during the PPN 
phase) . These studies did not address the question as to how 
the older, AGB spherical outflows would respond to "fresh", 
coUimated outflows from the PN central engine. 

In this paper we focus on a theoretical exploration of 
the role that changing mass-loss geometries from the central 
engine can play in driving changes in nebular morphology. If 
PPN drive coUimated outflows off the central engine but PN 
drive primarily spherical outflows then how will the nebula 
respond as it expands? The aim is to understand how neb- 
ulae change form as their hydrodynamic driving changes. 
Such process is of interest because it might help to under- 
stand the history of the central engine which may be a single 
star, a recently processed common envelope system or a bi- 
nary with an accretion disk. Thus there is the hope that one 
might be able to use morphological clues in an individual 
nebula to recover something of the history of the central 
object in terms of its mass-loss and the conditions which 
govern it. 

We note that there is a broader observational moti- 
vation for this study. The classiflcations of PPN and PN 
by morphological type appear to show signiflcant statistical 
differences. The general trends indicate that (i) around 70 
to 90% of PPN and young PN are b i polar or more complex 
l|Sahai fc Traugeilll998l : [Sahai et al.ll201ll l. (ii) about 20% 
of mature PN ar e round, 70% are elliptical and about 10% 
are bipolar (see iMiszalski et al.l l2009al : iParker et all I2OO6I : 
iLagadec et al.l I2OIII : the Planetary Nebula Image Catalog 



of Bruce Baliclo and references therein). There are how- 
ever considerable uncertainties in these classiflcations, e.g. 
(i) telescope resolution limitations; (ii) stellar population se- 
lection effects; (iii) differences in nebular class definitions 
(for example ring-like nebulae that are actually toroids in- 
cluded in spherical morphological classes); (iv) the well- 
established correlation between PN morphology and the pro- 
genitors' mass, i.e., bipolar and ring-like PN generally found 
at lower galactic sc ale heights, characteristic of higher-mass 
progenito rs (see e.g. Corradi fc Sch warz 19951: iKastner et al.l 
1I99& Manchado et al.ll201ll ). These factors, along with the 
fact that most attempts to classify PN by morphologi- 
cal type do not consider the ages of th e nebulae (but see 
Sah ai fc Traugeilll998l : ISahai et al.ll201ll for specific studies 
on young PN), complicate the interpretation of temporal 
changes in nebular morphology. 

If the imbalance between PPN and PN statistics holds 
up under further scrutiny then it implies that some objects 
may undergo a shape transformation as they expand. In par- 
ticular, one can imagine a scenario in which bipolar PPN 
hydrodynamically transform into elliptical or rounder PN. 
Such a transformation could be caused by different nebu- 
lar driving mechanisms during the early and late nebular 
evolutionary phases. 

Irrespective of issues of observational statistics, in this 
paper we take on the issue of morphological changes in 
PPN/PN from a theoretical perspective. Starting from a 
set of basic physical ingredients and using simulations, we 
aim to understand how shape transitions might occur dur- 
ing the evolution of an individual post- AGB system and how 
such morphological transitions depend on the parameters of 
the mass-loss process of the central star. The morphological 
evolution of single sources should be studied from first prin- 
ciples to better inform the interpretation of observations. 
Thus we consider the infiuence of an initial coUimated jet 
phase during the PPN and its effect on the previously de- 
posited circumstellar AGB wind. The jet phase is followed 
by a classical fast wind with a spherical geometry. Different 
distributions of the AGB wind's velocity and density are 
explored as are other parameters for the jet/wind interac- 
tions. Our choices for initial and boundary conditions are 
based on a reasonable balance between observational expec- 
tations and computational expediency. We note lSokerl (|l990|) 
has studied PN ansae using a similar wind/jet sequence, but 
here we carry out a more extensive study and explore a range 
of parameters across both PPN and PN evolution. 

This paper is organized as follows: in section [2] we de- 
scribe the methodology and numerical code used in the 
study, as well as our model and implementation of the AGB 
wind and post- AGB stellar outflows. The results of the sim- 
ulations are presented in section O In section 2] we compare 
our results with the general trends found by observations 
by means of synthetic emission maps and discuss how the 
results bear on issues of PN shaping. Conclusions are pre- 
sented in section [Sj 
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2 MODEL AND METHODOLOGY 

We carry out a set of Eulerian-grid numerical simulations to 
follow the interaction of coUimated jets and winds. The ini- 
tial and inflow conditions are meant to model nebular evolu- 
tion from the AGB to the mature PN phases. Wind dynam- 
ics are modeled and followed in 2.5 dimensions (cylindrical 
symmetry) using the equations of radiation hydrodynam- 
ics. The effects of optically thin cooling have be e n inc luded 
using the cooling tables of lPalgarno fc McCravl ([1973). We 
note that radiative cooling will play an important role in the 
shocks driven by the jet and also in the "rim" of swept up 
AGB material driven by the fast wind. Because of its high 
velocity the shocked fast wind should form a high tempera- 
ture "hot bubble" , Thb ~ lO'^ K, that will not cool effectively 
during the evolution of the PN. Observations have estab- 
lished, however, that only certain PN show such hot bubble 
emission and that the temperature is usually lower, by fac- 
tors of order 1 to 100, than that predicted by the jump 
conditions (e.g. iKastner et al.ll2008l . and references therein). 

The hydrodynamic equations are solved with the 
adaptive mesh refinement (AMR) numerical code As- 
troBEARl.Qj. In particular the Euler equations with cool- 
ing source terms are solved using a second-order MUSCL- 
Hanco ck shock capturing scheme and Marquina fiux func- 
tions (jCunningham et all [2003). While AstroBEARl.O is 
able to solve the equations of magnetohydrodynamics 
(MHD) and to compute several microphysical processes, 
such as gas self-gravity and heat conduction, we do not con- 
sider these processes in the present study. We justify these 
omissions at present on the basis that jets are likely formed 
via MHD processes at the central engine, but magnetic fields 
are no t likely to dominate global morphologies at nebular 
scales JHartigan et al.ll2007l '). 

The numerical domain of the simulations is 1 pc on each 
side. We use cylindrical coordinates (r, z) in a box that 
includes explicit calculations of the pole to avoid numeri- 
cal issues therqj. Thus our simulation extends from —rmax 
to Tmax, where Vmax = 0.5pc. Extrapolation, or outfiow, 
boundary conditions are used at r = ± 0.5 pc (above and 
below the z axis) and at z —1 pc. Reflective boundary con- 
ditions are imposed at the "equator" of the system (z = 0). 
We consider cells with radii R < r^, where r^ ~6000AU, 
to be the "nozzle" of the jets/wind and we apply inflow 
boundary conditions there (see below). 

A coarse grid of 128x128 cells is employed along with 
two levels of AMR refinement, each increasing the resolu- 
tion by a factor of 2. Thus the simulations attain a maxi- 
mum resolution of 512x512 with a AXmin ~403 AU. Typi- 
cal simulation flow times are of order 13000 yr which is about 
ten times longer t han previous s imulations of PN evolution 
fe.g. compare withlLee fc Sahail l2003. 2004: iLee et"ai]l2009l : 
lAkashi fc Sokeijl2008l ). We use BlueHiveFI. an IBM parallel 



^ Ihttp: //bearclaw. pas . rochester.edu/trac/astrobear/wiki 
^ We moved the bottom grid boundary away form the nebu- 
lar pole. Cells along this axis were evolved as grid cells, not as 
boundary (ghost) ones, by solving the equations of hydrodynam- 
ics. Hence we do not introduce any known numerical artifacts 
related to reflective boundaries on those cells. 
■* https://www.rochester.edu/its/web/wiki/crc/index.php/ 
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cluster maintained at the Center for Integrated Research 
Computing of the University of Rochester, to run simula- 
tions for about 20 hrs, using 16 processors. 



2.1 Outflow Phases 

Three outfiow episodes are considered in the simulations. 
The first on e, the "AGB wind" , con stitutes AGB initial con- 
ditions fe.g. lKnapp fc Morrislll985l ) set on the entire compu- 
tational grid. These conditions are such that the edge of the 
AGB envelope is located outside our computational domain. 
Our AGB wind/circumstellar environment has an isotropic 
radial velocity field of vagb =10kms~^ and a mass-loss 
rate of Magb = 1x10""' Mq yr~^. These values are repre- 
sentative of observed AGB wind properties (Hrivnak et al. 
1989; Bujarrabal et al. 2001). The AGB wind is set up with a 
uniform temperature of 500 K which best represents the in- 
ner region of the envelope. We would actually expect a tem- 
perature gradient, with the outer temperature being colder, 
however this simplification does not affect the current model. 

The second outflow phase, the "jet", begins immedi- 
ately upon the initiation of the simulation. A coUimated 
jet, with an opening angle of zero degrees, is injected 
at cells with radii < r^, (~6000AU). The jet's velocity, 
temperature and mass- loss rate are 200kms~^, 500 K and 
Mj =5xlO"^M0yr"^ (= 0.5 Magb), respectively. The jet 
is active for approximately 100 yr which is about 1% of the 
average total simulation run time. The launch and coUima- 
tion of our jets are assumed to occur in the "central erigine" 
(likely a binary or strongly magnetized stellar wind^) lo- 
cated at sub-resolution scales of order 10 AU. The jet injec- 
tion timescale is representative of expansion times (which 
are of order 500 yr) of some PPN and PN whic h can be 
directly imaged (e.g. see the PPN jet sample in iHuggins I 
l2007i . and references therein). Jet mass-loss rates are not 
well characterized so we have used an intermediate value 
(yet higher than those used in previous simulations, e.g. see 
iLee et al.l |2009| . and references therein) between the AGB 
and CSPN evolutionary stages, consisten t with both mo- 
ment um excesses observed in PPN flows (jBuiarrabal et al.l 
I2OOH ) and the observations of A 63 carried out by Mitchell 
et al. (2007). 

The third outfiow phase we con sider, the "fast win d" , is 
a standard central stellar fast wind (|Kwok et al.ll 19781) with 
a spherical velocity field continuously injected into the com- 
putational domain through the nozzle. In our three-phase 
models the spherical wind remains active throughout the 
remainder of the simulation after the jet is turned off. The 
fast wind has a mass-loss rate that decreases in time over 
~ 10* y r from 5xl0~^ to 5xl 0~^ Mq yr"\ following the re- 
sults of lPerinotto et al.l (|2004l . and references therein) . Once 
M reaches 5xl0~® Mq yr"'^ it remains constant at that rate 
until the end of the simulation. Simultaneously, the fast wind 
increases in speed from 200 to 2000 kms~^ in a manner that 
conserves wind ram pressure. This set of final fast wind con- 
ditions was chosen for computational expediency. Gas tem- 
perature in the wind at the nozzle is set such that the wind 



■^ We do not delve here on the actual physical causes of jets 
and other phenomena. We merely study the effect of sequences of 
plausible events on the PN morphology. 
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maintains a constant Mach number of 20 (again, for compu- 
tational expediency; the model results are not sensitive to 
this). 

We now review the suite of simulations carried out for 
the study. To isolate the effect of the jet we ran a series of 
two-phase simulations in which only the jet is launched into 
the AGB wind (Models 1 and 4, see Table 1). After the jet 
episode ends in these models no further gas is injected into 
the grid. Also, the simulations represented by Models 2, 2*, 
2b, 2b*, 5 and 7 are two-phase models in which the fast wind 
is injected directly into the AGB wind. The simulations rep- 
resented by Models 3, 3* and 6 are three-episode scenarios 
in which the fast wind is injected immediately after the jet, 
at time t ^ 100 yr (see Table 1). The simulation represented 
by Model 8 is almost the same as Model 3, except that the 
former has a wind-quiescent interlude that lasts ~ 400 yr 
between the jet and the fast wind episodes. The interlude 
represents an arbitrary timescale of order 4 times the jet's 
duty cycle (or active time) chosen to be long enough to al- 
low the jet driven cavity to evolve but short enough to be 
consistent with PPN/CSPN transition timescales. 

All simulations and their parameters, including the ad- 
ditional models that we present below, are summarized in 
Table 1. Evolutionary profiles of the outflow conditions are 
presented in Figure [T] We note again that the initial con- 
ditions are those of the AGB wind and set throughout the 
computational domain, whereas the conditions of both jet 
and fast wind episodes are set only on the gas located inside 
the nozzle {R < r^,, see above). 



2.2 Ionization front 

As the central star evolves from the AGB to the proto- 
white dwarf stages its surface temperature dramatically and 
rapidly increases and so docs its flux of ionizing UV photons 
jBalickl 1 19871 : iHoUenbach fc Tielens I Il997l : iPerinotto et al.l 
[2OOJ). Thus an ionization front will be driven away from 
the star into the surrounding circumstellar medium. Upon 
ionization the former neutral AGB envelope will be heated 
from temperatures of order T„eb ~ 10 A' to T^eb ~ 10^ K. 

The numerical description of ionization fronts with 
AMR requires the implementation of sophisticated meth- 
ods for radiative transfer. Efforts are currently underway to 
incorporate such methods into AstroBEAR2.0. In the cur- 
rent study we have implemented a simplified approximation 
which accounts for the thermal effect of ionization. When 
the fast stellar wind reaches a velocity of 2000 kms~^ we in- 
crease the temperature of cold gas (<500 K) to 10000 K. This 
temperature increase is imposed instantaneously across the 
entire computational domain. In this way we assume that 
the nebula is density bounde d and that the ionization front 
is R-typqj (see lGiulianilll98ll . and references therein), sweep- 
ing across the nebular domain in a timescale ti ^ R/c that 
is very short compared with the hydrodynamical evolution 
timescale th ~ R/v^ (where R, c and v^ are the radius of 
the nebula, the speed of light and the speed of the fast wind, 
respectively). In reality, however, ionization will not be as 
effective in the equatorial plane as towards the poles. 



^ In a R-type ionization front gas is adverted at supersonic veloc- 
ities and exhausts at a slightly lower, but still supersonic, velocity. 



Simulations that include this parametrization of ioniza- 
tion front passage are marked with an asterisk. For example, 
other than including the ionization-generated temperature 
increase, Model 1* is identical to Model 1, Model 2* is iden- 
tical to Model 2, etc. (see Table [1} . While this method is 
crude it allows us to test the potential of the pressure gen- 
erated via ionization to change the evolution of the overall 
nebular morphology at late states in its evolution. 



2.3 Additional models 

For completeness we have also chosen to carry out a se- 
ries of GISW models in which bipolar nebular morpholo- 
gies are expected to form in the PN phase (rather than 
through jets in the PPN phase). The GISW model re- 
lies on pole to equator density contrasts in an aspherical, 
toroidally shaped AGB envelope. Since toroidal AGB den- 
sity distributi ons have been adopted in GISW models be- 
fore (e.g. see iFrank fc Mellemal 119941 ) we have decided to 
run a series of models that include such density gradients as 
initial AGB wind conditions. Specifically, the AGB (ambi- 
ent/circumstellar) density distribution in Models 4, 5 and 6 



is the AGB aspherical density distribution used bvllcke et al 



( 1992, but see also references therein) and lFrank fc Mellema 
(1994 ), where 



Ptorus 



PAGB 

T(ff) 



(t)' 



F[e) = i-Q 



„/3cos(26))-/3-l 
R-2/9 - 1 



(1) 
(2) 



and PAGB = MagbI {^--kvagbt^)- 

This functional form produces a well characterized pole 
to equator density contrast. For example with a = 0.5 and 
/3 = 1 the above expressions introduce a pole to equator 
density contrast of 0.5 in the AGB wind. We have chosen 
these values in our study as they have been shown to produce 
butterfly type bipolar morphologies in pure GISW models 
(Frank & Mellema 1994). 

Model 4 is a two phase simulation and follows the inter- 
action only of the jet with the toroidal AGB wind. Model 5 
is also a two phase simulation but it follows only the fast 
wind interacting with the toroidal AGB wind. Model 6 is a 
three phase simulation tracking the jet followed by the fast 
wind interacting with the (preexisting) toroidal AGB wind 
(see Table [T|). 

Finally, we explore the role of velocity gradients in the 
AGB wind. We have carried out a simulation in Model 7 in 
which Vagb = V{6) = /(r, z) * Vo where the constant radial 
velocity field Vo is multiplied by the function 



fir,z) 



1 



^-tan-U\r/^\)V0.3^ 



(3) 



which results in a latitudinal decreasing expansion velocity, 
with a polar-to-equatorial contrast ratio of 2. A continuous 
fast wind is then driven into this aspherical AGB wind. As 
we shall see in the next section such a velocity gradient does 
not greatly affect the morphology of the resulting nebula. 
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Table 1. Simulations and parameters. 



Simulation 


AGB wind 


Jet 


Fast wind 


Fast wind 


Ionization 


name 


distribution 


duration 


duration 


max. speed 


front 






[xlOOyr] 


[xlOOOyr] 


[xlOOOkms-i] 


passage 


Model 1 


spherical 


1 


0.0 




n 


Model 1* 


spherical 


1 


0.0 




y 


Model 2 


spherical 





13.0 


2 


n 


Model 2* 


spherical 





13.0 


2 


y 


Model 2b 


spherical 





16.3 


1 


n 


Model 2b* 


spherical 





16.3 


1 


y 


Model 3 


spherical 


1 


10.7 


2 


n 


Model 3* 


spherical 


1 


10.7 


2 


y 


Model 4 


toroidal p 


1 


0.0 




n 


Model 5 


toroidal p 





3.8 


2 


n 


Model 6 


toroidal p 


1 


6.0°^ 


2 


n 


Model 7 


aspherical v 





13.0 


2 


n 


Model 8^ 


spherical 


1 


18.0 


2 


n 



" This time is shorter than the one of Model 3 (10.7 kyr) because the toroidal AGB 
wind in this model funnels the nebula's lobe and thus it reaches the computational 
domain's boundary faster than the lobe in Model 3. 

Model 8 includes a relaxation interlude of ^^ 400 yr between the jet and the wind 
episodes, and a computational domain of 2pc^, instead of Ipc^ as in other simula- 
tions. 



3 RESULTS 

The results of our simulations are presented via logarithmic 
false color gray scale maps of the gas density with overlaid 
velocity vectors in Figures [JHH Panels in these figures are 
arranged such that columns correspond to different simula- 
tions and rows show time increasing from top to bottom. 



3.1 Structure within spherical AGB winds 

Jet Models: In figure[2] left column, we show the results of 
Model 1, in which a short coUimated jet phase is driven into a 
spherical AGB wind envelope. The figure shows the jet driv- 
ing a bow shock which bounds a narrow-waisted cavity. The 
cavity wall is defined by a thin dense "rim" and its bipolar 
geometry persists for the entire duration of the simulation. 
The lateral width of the lobe is the result of the thermal pres- 
sure within the cavity as shocked material is forced laterally 
out of the region of the jet head. Thus the lobe's width-to- 
length ratio is determined by the sound speed in the AGB 
shell and the speed of the jet head. We note that the cool- 
ing time of the post-shock gas, t c '-^ T/(nA(T)), wh ere A(r) 
is the cooling function (Dalgar no fc McCravHigTa ). is short 
relative to the hydrodynamic time (which is the ratio of the 
computational domain length over the gas sound speed). 
Thus all post-shock flows in this simulation are relatively 
isothermal. 

After about 3500 yr, the cavity lobe develops a narrow 
waist. This timescale is of order the expansion time of ma- 
ture, but not old, PN. The pinch at the waist is the result 
of the higher density (and, therefore, higher ram pressure) 
of the AGB shell at small radii. In section 14.11 we analyze 
the evolution of the aspect ratio of some of our model neb- 
ulae. Here we note the aspect ratio of the cavity's bounding 
bipolar lobes increase monotonically in time. 



Fast Wind Models: In figure^ right column, we show 
the results of Model 2, in which a "classical" fast wind is 
driven into a spherical AGB wind envelope. The growth of 
the nebula follows the analyt ical predictions of Kwok et al. 
(1978, but see also lStute fc Sa hai 2006) and agrees well with 
Schonberner, Jacob, fc Steffen (2005), whose 1-D models in- 
clude much more detailed physics that we could include in 
our multidimensional models. The slow dense AGB wind is 
quickly overtaken by the spherical fast wind which is ini- 
tially 100 times sparser and 20 times faster. An inward fac- 
ing shock decelerates fast wind gas which is then heated to 
temperatures of T ~ 10^ — lO^K; a spherical hot bubble is 
formed. Such high temperatures have never been observed 
in PN hot bubbles, and, hence, a number of mechanisms 
ha ve been proposed to moderate their te mperatures (e.g . 
see lSoker et al.ll2010l . and references therein: iLi et al.ll201ll ). 
The thermal pressure of the hot bubble acts as the piston 
driving a shock into the AGB envelope. AGB gas swept up 
by this shock cools rapidly and is compressed into a dense 
rim expanding with a velocity of Vrim ~ 20kms^^. In this 
way an isotropically expanding PN is formed, just asde- 
scribed by the simple analytical model of lKwok et al.l ([1973) ■ 
The nebula reaches a radius of '^ 0.4pc after expanding for 
about 5000 yr. 

We note that observations (Bruce Balick 2011, private 
communication) often find that AGB shells expand with an 
outward radial gradient (i.e. vagb oc Ragb)- This is not 
the case in our AGB initial conditions (section 12. ip . How- 
ever, we carried out additional test simulations (not shown) 
and have found that the morphology results of the present 
paper do not sensitively depend on the functional form, 
VAGB ~ constant vs vagb oc Ragb , of AGB environments. 
This happens, in part, because the velocity of the jet and 
fast wind is faster than vagb by at least one order of mag- 
nitude. 
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Figure 1. Outflow evolutionary profiles of Models 1-8. The initial conditions (t^O.OOl yr) are those of the AGB wind and arc set 
throughout the computational domain, whereas the conditions of both jet and fast wind phases are set only on the gas located inside 
the nozzle. See Table 1 for details. 



Jet and Fast Wind Models: In figure O left column, 
we show the results of Model 3, in which the coUimated jet 
phase is followed by the spherical fast wind. Initially, the 
collimated jet forms a bipolar cocoon at the center of the 
AGB envelope as in Model 1. Once formed, the jet-driven 
cavity is then shaped from within by the isotropic fast wind. 
The fast-wind shock quickly forms, filling the initially bipo- 
lar cavity with high temperature gas. The contact surface 
of the jet-driven cavity is pushed against the slow denser 
AGB envelope by the internal pressure of the fast wind. Be- 
cause the jet has already shaped the AGB wind, an elliptical, 
rather than a spherical, dense rim is eventually formed. 

The bubble attains an aspect ratio of ~ 2 which is 



achieved early in the evolution and remains roughly constant 
as the nebula expands. This simulation demonstrates one of 
the principal points of our study: An initially bipolar out- 
flow driven by a strongly collimated jet can be transformed 
into an elliptical nebula at longer times and larger scales by 
a subsequent isotropic, fast wind. Thus nebular morphology 
is not a given of initial conditions but can change over time 
as the wind driving conditions change. 

In most of our models the transition between jet and 
wind occurs instantly. In Model 8 (figure [5l right column) 
however, we show the results of a simulation in which we 
have modified the conditions in Model 3 to include a quies- 
cent period between the jet and the wind. In the evolution of 
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real PN this gap could represent a phase during which energy 
injection from the CSPN, or binary system, may be absent. 
As discussed in the previous sections time-scale estimates 
of such episodes are of order the jet duty cycle, t ~ 100 yr 
(see section [2TT]) . The simulation shows that, once again, the 
jet phase leads to the formation of a cavity that is initially 
bipolar. This fossil structure then expands quasi-ballistically 
for ~400yr, reaching scales of order 0.05 pc. The fast wind 
then begins at time t = 500 yr, and the evolution of the in- 
teracting winds is essentially similar to that of Model 3. 

We note that in this simulation the rim of swept up 
AGB material maintains an aspect ratio close to 2.3 which 
means that once the fast wind begins the nebula expands 
homologously (see Figure [6] and section [4^ . One difference 
between this simulation and Model 3 is the late time for- 
mation of a dense knot along the long axis of the elliptical 
rim. The number density and flow speed in the knot are of 
order 1500 cm~^ and 30kms~^, resp ectively, which are not 
unlike those in some FLIERS (see e.g. lBalick et al.lll998l . and 
references therein). The formation of the knots appears due 
to the radiative collapse of jet material as it interacts with 
the AGB wind. In contrast, when the fast wind immediately 
follows the jet phase (Model 3) the flux of mechanical energy 
into the polar regions flattens the jet material and keeps it 
from forming a ballistically expanding dense knot. 



3.2 Structure formation in non-spherical AGB 
shells 

The nebular dynamics changes significantly when the AGB 
envelope begins with a toroidal density distribution. The 
results of our studies in this context are summarized in Fig- 
ure[4l We note that in these models the jet propagates along 
the symmetry axis of the toroidal AGB envelope. 

Jet Model: In figure [3l right column, we show the 
results of Model 4, in which the jet alone drives into an as- 
pherical AGB wind. Once again the jet drives a fossil bipo- 
lar cavity which develops a narrow waist at ~ 3000 yr com- 
pared to 5000 yr for the spherical AGB sheU (of Model 1). 
The lateral expansion of the lobe is, once again, stronger in 
lower density regions where the ram pressure of the initial 
AGB envelope is easier to displace. The lobe that develops 
in this case has a narrower waist than in models with no 
pole to equator density contrast. The dynamics of the evo- 
lution is, however, not significantly different from the results 
of Model 1. This means it would be difficult to determine 
from the morphology of a PPN whether or not the progen- 
itor AGB wind had a significant pole to equator density 
contrast. Jets, if present, dominate the morphology only if 
there is no fast wind or if the fast wind has less momentum. 

Fast Wind Model: In figure IH left column, we show 
the results of Model 5, in which only a fast light isotropic 
wind interacts with the dense slow toroidal AGB envelope. 
This is the classic GISW model that has been explored be- 
fore (e.g. Icke et al. 1992; Frank & Mellema 1994). The Fig- 
ure shows that, as expected, the swept-up AGB shell be- 
comes elliptical at early times and then becomes more sig- 
nificantly bipolar as the hot bubble drives the shock to larger 
radii in the AGB wind along the polar axis. At late times the 
dense rim takes on a butterfiy morphology as would be ex- 
pected for these initial conditions. Thus this simulation and 
that of Model 2 (spherical AGB wind) recover the results 



of previous numerical experiments in the GISW and ISW 
models (e.g. iFrank fc Mellemal 1 19941 ) and thereby confirm 
the consistency of our numerical simulations. 

Jet and Fast Wind Model: In figure [4l right column, 
we show the results of Model 6, in which the jet is followed 
by the fast isotropic wind as they interact with the toroidal 
AGB envelope. Once again an initially bipolar cavity and 
rim are created by the jet. Once the fast wind begins how- 
ever the evolution differs significantly from Model 3. Rather 
than changing the morphology from bipolar to elliptical, the 
action of the fast wind in this case is different; the evolution 
is dominated not by reshaping the jet driven cavity but by 
interacting with the toroidal AGB wind. The long term evo- 
lution of this model leads to a butterfly shaped rim which is 
quite similar to that seen in the previous simulation in which 
no jet was included. Thus it appears that the presence of a 
preexisting density contrast in the AGB wind is more impor- 
tant for the evolution of morphology than the presence of a 
short-lived jet phase as no change from bipolar morphology 
to elliptical morphology occurs. We note a l arge fraction of 
optical PPN have dense dust waists (see e.g. iHuggins 1120071 . 
and references therein). Hence the above conclusion would 
suggest that a large fraction of mature PN should also be 
bipolar. Since this is not the case it is possible that such 
dust shells may not be extended enough to influence PN 
morphology once the fast wind begins. 

Asymmetric AGB Velocity: In figure[S] left column, 
we present Model 7, in which the fast wind expands into an 
AGB envelope with an aspherical velocity (rather than den- 
sity) distribution. In our initial conditions the AGB wind 
expands two times faster toward the pole than towards the 
equator, following equation ((3]). The hot bubble responds 
to the reduced AGB ram pressure Prp{6) — pagbVagb{6)^ 
along the equator leading to high shock velocities at these 
latitudes. As the figure shows the resulting morphology is el- 
liptical. This simulation demonstrates that if the early evo- 
lutionary phase of the PN produce a r apidly expanding equa - 
torial flow rather than a bipolar jet (JNordhaus et al.ll2007l l. 
then an elliptical nebula might be the long term result. 



3.3 Effect of Ionization Induced Temperature 
Increase 

We note finally that we have compared the evolution of those 
models with our simple estimate of ionization against the 
ones with no change in temperature (see Table 1, last col- 
umn). We find only small changes, of order 5%, in the dy- 
namical behavior of simulations with a temperature increase 
(from T < 500 K to T = lO'' K). 

Our models are relevant to cases in which the ioniza- 
tion front remains R-type (see section 12. 2p until it leaves 
the AGB material. We are not able to see details of the 
ionization front's passage on the gas in terms of driving in- 
stabilities which can fragment the dense rim when the fronts 
stall in the nebular gas. We note that even 1-D models show 
that the trapping of th e ionization can have a n effect on 
the evolution of the rim (jSchonberner et aLlbOOST ). These ef- 
fects should not however affect our conclusions about global 
changes in morphology. For AGB envelopes that begin with 
little pole to equator density contrast the trapping of the 
D-front will occur at all latitudes and its changes will affect 
the entire rim. Thus ionization effects should not change 
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Figure 2. Dynamical evolution of the gas for different simulations that started with a spherical AGB wind distribution. These are 
logarithmic false colour grayscale maps of the number density. Columns correspond to different simulations and time increases downwards. 
Arrows indicate velocity. See Table 1 for details on the model parameters. 



our conclusions about the morphological shift from bipolar 
to elliptical nebulae as they will not change the nature of hot 
bubble pressure gradients driving the rim from the inside. 
We leave a detailed study on this matter for future work. 



4 DISCUSSION 

The main results from our simulations are the following. A 
young bipolar PPN transforms into an older elliptical PN 
when an initial spherical AGB envelope interacts with a 
short duration jet and then with a spherical post-AGB fast 
wind. Once mature PN become elliptical they do not even- 
tually change into more spherical nebulae over relevant time 



scales. If the initial AGB envelope is toroidal then butterfly 
nebula occur with or without a PPN jet phase. 



4.1 Aspect ratio evolution 

We follow the morphological transition of our model nebulae 
formed from spherical AGB envelopes (Figure [21 Figure [3]- 
left and Figure [5]-right) using plots of the rims' aspect ratio, 
e(t) = Lma{t)/Lmi{t), as a function of time. This is the 
major-to-minor axis ratio of a rectangle that would contain 
the principal bright structure of the nebulae. The aspect 
ratio evolution profiles are shown in Figure [B] It is clear that 
their gradients are related to the outflow phase sequence, or 
history, in the simulations. 



© 2012 RAS, MNRAS OOO.ITHTSi 



Simulating the Morphological Evolution of PN 9 



"f-T -1.1 -<f-r fl-# I, J 



1,4 S,9 



MOs. 1 yr 



0.4- 



#.i- 



O.0-, 



iftffff/ff'f'**'^ f ^ * ^ * ^ i» -f,^ * * * ^ ^ n' ^ * 




^u;.; yi 



1 1 J /# ' * • 

i > f i f t ' ,',' 

ltttf*/r 
y 1 1 f f * *f 
\ I r r 4 * * r* Ji»*^'* 

' r f f ' • .» JTrr^ _ ^ - - 
^ >-'- - ■ _ — — 




T^O.^yr 



^"'5.2 4 ■ 



I f t f f f ft ^ ^ ,' f f r^J' J*J'J> .-.jr,rj'^^,r^jr ^^ ^ ^ jr 

t ft! ! ' f Jt * f * > 




Mattel St 9phsritBl AGS vniui ~i joJ -) /lui iL'ffu^ 



0.2 0.4 0.6 

pc 



o.s 

Madd 4; tamiiai AGB * jet 



Figure 3. Dynamical evolution of the gas. These are logarithmic false colour grayscale maps of the number density. Columns correspond 
to different simulations and time increases downwards. Arrows indicate velocity. See Table 1 for details on the model parameters. 



The spherical PN rim that we see in Model 2 has a 
constant e(i) of unity (thick solid line). The bipolar rim in 
Model 1 has an e(i) that increases monotonically as the ob- 
ject expands (thin dashed line, Figure O. The expansion of 
this lobe is dominated by the momentum of the jet during 
the early evolutionary phases, hence we see a steep increase 
in e(t). Note that the slope decreases in time because the jet 
turns off. If the jet were active we would see a constant slope. 
The aspect ratio of the lobes is (of course) also proportional 
to the duration of the jet injection phase. 

In contrast, the rims in Models 3 and 8 (the thin solid 
and thick dashed lines in Figure |6] respectively) show quite 
different evolution of e{t). For these simulation e(i) increases 
steeply during the early evolutionary phases in which the jet 
dominates the object's expansion. This corresponds to the 



PPN phase. Any further increase of the aspect ratio is then 
quickly suppressed by the effects of the spherical fast wind. 
The profiles reach values of about 2 and vary only modestly 
afterwards. This means the expansion of the rim becomes 
homologous (section l4.2p . This is the most important point 
of the work. 

The thick dashed line corresponds to Model 8 and shows 
a short increase/decrease in aspect ratio that occurs be- 
tween the early evolution and the later phases when the 
aspect ratio becomes constant. This feature occurs during 
the quiescent evolutionary phase, i.e. the PPN-PN transi- 
tion. Comparing with the aspect ratio profiles of Models 1 
and 3 we see that the height of the bump in Model 8 is 
proportional to the duration of the quiescent interlude (see 
end of section 12.11) . An interesting question here is that of 
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Figure 4. Dynamical evolution of the gas for different simulations that started with a toroidal AGB wind density distribution. These are 
logarithmic false colour grayscale maps of the number density. Columns correspond to different simulations and time increases downwards. 
Arrows indicate velocity. See Table 1 for details on the model parameters. 



the dependence of e{t) on the duration of this interlude. To 
answer this, additional simulations (designed to explore a 
set of values for the duration of the interlude) would have 
to be carried out; we leave this for future research. 

Our calculations therefore demonstrate how the mor- 
phology of PPN and PN is correlated with the coUimation 
parameters of the wind from the central object. The pres- 
ence of a brief, hypersonic and heavy jet prior to the interac- 
tion between the AGB envelope and the spherical fast wind 
changes the shape of the resulting nebulae. The effects are 
important during both the early and late phases of nebular 
expansion. Bipolar PPN are formed and then turn into el- 
hptical PN, in broad agreement with the general statistical 
trends of the observed morphologies of PN. 



4.2 Velocity field 

We show the velocity field of gas in our simulations in Fig- 
ures 2, 3, 4 and 5. Vector density, location and length are 
constant, but vector inclination and grey scale are related 
to the local direction and speed of the flow, respectively. 
Arrow properties were chosen in such a way as to stress ve- 
locity differences between the polar and equatorial parts of 
our model nebular rims. A detailed study of the velocity 
field distribution of material inside the lobes, as well as in 
butterfly-shaped rims, is beyond the scope of this paper. 

The following findings are consistent with the nebular 
aspect ratio evolution profiles of Figure [6] The left column 
of Figure [2] shows that the tip of the bipolar rim formed 
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Figure 5. Dynamical evolution of the gas. These are logarithmic false colour grayscale maps of the number density. Columns correspond 
to different simulations and time increases downwards. Arrows indicate velocity. See Table 1 for details on the model parameters. 



in Model 1 maintains a longitudinal (horizontal) velocity, 
Vz, which is close to 75knis~^. Lobe material located be- 
hind the tip shows a v^ increase proportional to the distance 
from the star, z. This is consistent with the "Hubble flow" 
kinematics observed in some PPN outflows (Balick & Frank 
2002; Bujarrabal et al. 1998; Olofsson & Nyman 1999). In 
contrast, the velocity of the rim's waist (Figure [21 left) is 
mostly transverse, Vr, slower (darker) and quite similar to 
that of the AGB wind material located at larger radii. We 
see a similar velocity distribution at the tip and waist of 
the bipolar rim formed in Model 4 (Figure [31 right column). 
I.e. a fast (light-gray) tip with Vr,Vz ~ 0,100 kms~^, and a 



slower (dark-gray) waist with Vr,Vz ~20,0kms~^. In this 
case however most of the lobe material, not just the compo- 
nent behind the rim's tip, also shows a v^ increase propor- 
tional to the distance from the star. 

The left column of Figure [3] shows that the elliptical PN 
rim of Model 3 has a consistent, progressive speed increase 
(dark to light gray) of a factor ~ 2, form the waist to the 
tip along the rim. This velocity field goes from transversely 
dominated at small z to longitudinally dominated at large z. 
The object expands homologously. In this case we do not see 
any correlation between the position on the z-axis and the 
speed of the lobe material. We see similar velocity distribu- 



© 2012 RAS, MNRAS OOO.ITHTSl 



12 M. Huarte-E spinas a et al. 



g 

o 
a> 
a. 



Ml: AGB->jet 

M2: AGB->fastwind 

M3: AGB->jet->fast wind 

M8: AGB->jet->gap->fast wind 



5000 



10000 15000 
Time [yr] 



20000 25000 



Figure 6. Time evolution of the nebular aspect ratio, e(t). The 
1/— axis shows the major-to-minor axis ratio of the rims (the prin- 
cipal bright structures). Profile gradients correlate with the out- 
flow phase sequence, or history. 



tions in the elliptical PN rims of Models 7 and 8 (Figure [5]). 
The tip-to-waist speed ratio [vz{r — 0)/vr{z — 0)] of the 
rims in Models 7 and 8 is '--^ 1.5 and 2, respectively. Also, 
we consistently find velocity vortices inside the lobes of our 
model elliptical PN. I t is worth not i ng th e semi-analytical 
kinematic analysis of ISteffen et al.l (|201ll ) which suggests 
that in general PN rims asymptotically evolve towards a 
homologous expansion, in agreement with our findings. 



4.3 Synthetic observations 

The simulations described in section 3 yield detailed pre- 
dictions of the density structures of PN (as well as their 
detailed, if simplified, temperature structures) that, in prin- 
ciple, can be used to generate synthesized emission maps 
for comparison with the main observed PN morphologies 
across the electromagnetic spectrum. As a first step in this 
direction, we synthesized emission measure distribution im- 
ages (i.e. the integral of p^ along the line of sight) for 
Models 1, 2, 3 and 8 using the reconstruction tool Shape 
(|Steffen fc Lopej I2OO6I) a nd fo llowing the procedure de- 
scribed by ISteffen et al.l l|200g ). We choose these models 
because their shapes depend only on their outflow histo- 
ries. We use logarithmic grayscales and an inclination of 
90° between the polar axis of the nebulae and the line of 
sight (edge-on). Note that in constructing these maps, we 
have not attempted to account for spatial variations in gas 
temperature or ionization state; because much of the gas 
in these simulations is at time average mean temperatures 
of order 10^ K, the synthetic images in Figure 7 likely best 
approximate the appearances the nebulae would have in op- 
tical emission lines (e.g., reflected starlight, emission lines 
arising in shocks and recombination lines). A more detailed, 
multiwavelength treatment of the synthesized emission mor- 
phologies, which is beyond the scope of this paper, will be 
the subject of future work. 

With the above caveats in mind, these initial synthe- 
sized emission maps illustrate the clear relationship between 
the emission morphology of a PN and its CSPN outflow his- 
tory (i.e. coUimated jet or isotropic fast wind). A bipolar 



object is formed by the interaction of a dense slow AGB 
envelope with a brief phase of coUimated heavy jet ejection 
(Model 1, Figure 7). The interaction of a light isotropic fast 
wind with a spherical AGB envelope yields a spherical neb- 
ula (as in Model 2 as well). The PN rim in these cases 
shows an emission distribution with a smooth interior sur- 
face brightness and limb-brightened edges. 

When all three outflow phases are present, on the other 
hand, we see young (not-narrow-waisted) bipolar nebulae 
(Models 3 and 8, top row. Figure 7) transformed into older 
larger elliptical nebulae (bottom row). The process is par- 
ticularly evident in Model 8 and affirms the behavior seen 
in the aspect ratio profiles (Figure [G]). Such morphological 
transition is consistent with the general trends found by ob- 
servational studies. 

The elliptical rims in cases 3 and 8 have similar mor- 
phologies. However, their synthetic emission during the PN 
phase is significantly affected by the presence of the quies- 
cent episode simulated only in Model 8. Figure 7 (bottom) 
shows that the nebula in Model 3 has an emission distribu- 
tion with only modest variation in brightness along the rim. 
In contrast, the PN nebula in Model 8 shows two bright 
knots along its polar axis, located symmetrically with re- 
spect to the center. Examination of the simulation data 
shows a density contrast of order 100 between the main 
part of the nebula and the bright knots. It is noteworthy 
that these features bear some morphological resemblance to 
the FLIERS that are seen in some elliptical PN (see end of 
sectionl3.1l). 



4.4 



PPN/PN morphological changes and 
progenitor stars 



It has been suggested that the outflow mass-loss properties 
of PN, in particular their geometry, are affected by both 
the charact eristics and multiplicity of their progenitor stars 
( Ibcn I99I lLivid[T993: Soker 1997; Soker 1998, and refer- 
ences therein; Nordhaus fc Blackman 20061: Nordhaus et al.l 



I2OO7I : iDe Marcdl2009l : iDe Marco fc Sokeill201ll ). Where do 
our models fit in this picture? In other words, which binaries 
are expected to produce both spherical AGB envelopes and 
jets which could then lead to the morphological change we 
see in Models 3 and 8 (from bipolar to elliptical)? 

Although we have a reasonable understanding of how 
dust-driven winds can explain the high mass-loss rates ob- 
served during the AGB superwind phase (e.g. Morris 1987; 
Lagadec & Zijlstra 2008) questions remain open, in partic- 
ular regarding the geometry change that must take place 
during the AGB-to-PPN transition. Today, it appears un- 
likely that a single star may power the evolution of an AGB 
star that results in a non-spherical PN (e.g. Nordhaus et al. 
2007), although it is possible that in the future we may un- 
derstand how a dynamo could be sustained in an AGB star 
and generate a super-wind phase that departs from spher- 
ical symmetry (Blackman et al. 2001). The type of binary 
that could give rise to the simulated ejection phases can be 
surmised from basic physical arguments and simulations. 

Companions that are located within 2-3 AGB star radii 
can be captured into a common envelope (Nordhaus et al. 
2010). The subsequent in-spiral causes the AGB envelope to 
be ejected (likely in a toroidal geometry; Sandquist et al. 
1998) and the post-AGB star to emerge with a much re- 
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Figure 7. Synthetic emission maps of the model nebulae. The inclination angle o f the images is 90°. These are emission measure 
distribution images (i.e. the integral of p^ along the line of sight) produced using Shape l ISteflFen fc L6pe2l2006l) . Grayscales are logarithmic. 



duced radius and much higher temperature. Observed sys- 
tems (e.g., De Marco 2009) and simulations (e.g., Sandquist 
et al. 1998, Ricker and Taam 2008, Passy et al. 2011) in- 
dicate that the post-common envelope primary will be a 
star with a radius within a solar radius and may bypass the 
PPN phase altogether, emerging primarily as a bipolar PN 
(althou gh some post-common envelope PN are actually el- 
lipticah lMiszalski et al.ll2009al lbr). It is not clear whether a jet 
would ensue, although highly collimated structures are ob- 
served in post-common envelope PN (e.g., Abell 63; Mitchell 
et al. 2007). 

For some post-CE PN with jets, the jet is kinematically 
older than the main PN and may have started before the 
bulk of the envelope departed the system (e.g. A63, Mitchell 
et al. 2007) or ETHOS 1 (Miszalski et al. 2011). For these 
cases Models 4, 5 or 6 may be fitting. In at least one case, 
however, this is the other way around and the jet is younger 
than the PN (NGC6 778, Guerrero & Miranda 2012). 

A class o f bina ries exists that has larg er final sepa rations 
(100-500 Rrr,: lvan W inckcl 2003; van Win ckel et al.| [2009'). It 
is unknown if a common envelope resulted in their forma- 
tion. Such post-AGB primaries are almost never accompa- 
nied by a PPN. In a couple of instances bipolar outflows 
are seen (e.g., the Red Rectangle; Cohen et al. 2004, and 
references therein) . We note that circumbinary disks are ob- 
served around all of these system. These binaries may never 
develop a proper PN due to their envelope being replenished 
by fall-back of disk material, which could prevent a temper- 
ature increase of the central star. It is not clear if any of our 
models might mimic this evolutionary channel. 

Other binary types that could fit our models can be 
wider binaries, those with initial separations larger than a 
few stellar radii (5-10 AU). These can f ocus the AGB win d 
into an equatorially concentrated torus (JEdgar et al.ll2008h . 
accrete matter and blow jets either before or after AGB 
departure. Any of our models could simulate such scenario. 
If the companion is lower mass, both the focussing action 
and the jet blowing may be weak. The star might evolve as 
if single then. 

If a companion is captured tidally by an AGB star and 



it has considerable mass (larger than a few tenths of a solar 
mass), it will sy nchronise the AGB stellar rotation with the 
orbital motion (|Counselmanlll973r ). When this happens the 
AGB star is spun up and may eject an equatorially concen- 
trated wind. A common envelope would be avoided because 
the tidal capture would be slowed down. Jets may be blown 
for a period during or after the AGB envelope ejection. An 
example of this scenario could be Model 6. 

Finally, if a companion caught into a common envelope 
is low mass (such as a brown dwarf or even less massive) 
and tidally destroyed within the common envelope, it may 
form a disk around the core and blow jets. In this scenario 
a jet may precede a possibly toroidal AGB wind. 



5 CONCLUSIONS AND SUMMARY 

We have carried out a series of 2.5D hydro simulations de- 
signed to study how changing the mass-loss properties of 
PN progenitor stars affects the evolution of nebular mor- 
phology. The stellar mass-loss evolution has been followed 
starting at the AGB phase, evolving through the PPN phase 
and finishing at the mature PN phase. 

Our simulations show that a young bipolar PPN trans- 
forms into an older elliptical PN when an initially spherical 
AGB envelope interacts with a short duration jet, driven 
from the central star and followed by an isotropic fast wind. 

Thus we have demonstrated how changing mass-loss 
geometries from the central engine can drive evolutionary 
changes in nebular morphology. Our theoretical results com- 
plement GISW studies. A fast spherical CSPN wind driving 
into a slow spherical AGB wind has been thought to create 
a spherical nebula. We have shown that if this interaction is 
preceded by a collimated PPN jet phase then an elliptical 
PN will result. A fast spherical CSPN wind driving into a 
slow aspherical AGB wind has been thought to create an 
aspherical nebula and our results show that including an 
earlier PPN jet phase does not affect the eventual butterfly 
morphology that results. 

Our theoretical result that bipolar PPN can evolve into 
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elliptical PN is consistent with general observational trends 
suggesting that PPN appear to favor bipolar morphologies 
while mature PN appear to be more elliptical. This result 
has to be considered along side of the complementary obser- 
vational trend that bipolar PN are found at typically lower 
Galactic scale heights than ellipticals, suggesting a correla- 
tion betwee n bipolarity and higher-mass (shorter lifetime 
progenitors llCorradi fc Schwarj|l995l : iKastner et al.lll99i 
iManchado et al.lbOlll 'K 

We also have found that once mature PN become ellip- 
tical they do not evolve to become spherical nebulae over 
the relevant PN lifetim es, so the origin of mature spheri- 
cal P N (19% of all PN; iParker et al]|2006l : iMiszalski et"ai] 
120081 ) does not seem to follow the initially aspherical path- 
way studied herein. 

The aspect ratio evolution of our model nebulae sug- 
gests that bipolar PN with projected aspect ratios >4 
may result from CSPN which during the post-AGB phase 
develop a brief jet, but no fast spherical winds, or rela- 
tively weak ones (e.g., symb iotic Mira systems like R Aqr, 
Mz 3, Belczviiski et al.ll2000l . and references therein; Mz 3, 
ISchwarz et al.lll992l ). 

When the initial jet drives into a toriodal AGB density 
distribution the result is a bipolar PPN that evolves into a 
mature bipolar nebula with a narrow waist. The toroidal 
AGB density distribution in this case dominates the dy- 
namics and shape of any subsequent outflow phases. At late 
times the asymptotic geometry of this nebula is a traditional 
butterfly nebula when seen in projection onto the sky. We 
also flnd that the interaction of a spherical post-AGB fast 
wind with an AGB envelope which has an aspherical velocity 
field yields a homologously expanding elliptical nebula. 

One of our simulations included a quiescent interlude 
between the jet and the fast wind episodes. We find the cor- 
responding synthetic emission map shows small bright fea- 
tures which resemble observations of axial knots or FLIERS 
within PN. A future avenue for these studies would be to ex- 
plore when and how knots form as the result of jet material 
collecting at the head of the flow. 
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